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Abstract 

We present the crumpling transition in three-dimensional Euclidian space 
of dynamically triangulated random surfaces with edge extrinsic curvature and 
fixed topology of a sphere as well as simulations of a dynamically triangulated 
torus. We used longer runs than previous simulations and give new and more 
accurate estimates of critical exponents. Our data indicate a cusp singularity 
in the specific heat. The transition temperature, as well as the exponents are 
topology dependent. 



PACS: 05.90, 64.60A, 64.60F, 68.10, 87.20^, 12.90 



1 



The Crumpling Transition of Dynamically Triangulated Random Surfaces 



2 



1 Introduction 

In this paper we are concerned with the statistical mechanics of surfaces. A possible 
approach in space dimensions D > 1 is to discretize the surface using a triangulation 
[1-4]. For this we need to specify the Hamiltonian Ti. The surface S is replaced 
by a simplical triangulation T, specified by the number of nodes A^, of links and 
triangles, and the X-coordinate field by the coordinates X of the nodes. The metrical 
fluctuations of the manifold are modeled by summing over triangulations induced 
by link-flips [5-7]. The Hamiltonian is now choosen, such that the partition function 
is not dominated by configurations with spikes. In order to surpress these spikes one 
adds a term with extrinsic curvature. As a function of the extrinsic curvature the 
model shows a transition (crumpling transition) at finite rigidity of the surface [8-18]. 
In a previous paper, we discussed the extrinsic and intrinsic geometrical properties 
of dynamically triangulated random surfaces (for example the Hausdorff- dimension. 



spectral and spreading dimension) above, below and at the transition ||32 
The partition function of the above model can be written as 



N-l 

Zm= d'^xJ n d'^X.e-^ (1) 
1=1 



where the translational mode is integrated out. The Hamiltonian Ti is defined as 



N 

^ ^ ^ ^ 

Hg 

The Gaussian part of the Hamiltonian Tig is a sum over the positions X in embedding 
Euclidian space of all nearest neighbours nodes, i.e. all links of the triangulation. 
We shall use (3 = 1 because of the rescaling invariance. 

He is an edge extrinsic curvature term [19-26]. Z^a^.a^ denotes a summation over 
all adjacent triangles which share an edge and h/^- ■ h/\. is the scalar product of the 
vectors normal to a triangle. 

The third part of the Hamiltonian Tim is the discretization of the square root of the 
metric. (Tj denotes the number of nearest neighbours of node i. a depends on the 
measure. We used a = D/2 for the simulations in D = 3 dimensional embedding 
space. 

Here we investigate such a model numerically (Monte Carlo simulations), apply- 
ing finite-size scaling ideas. These ideas have proven successful in other contexts for 
the analysis of phase transitions. We expect that the above model shows a phase 
transition at some critical value of the coupling A. There are several questions which 
we want to address: What is the order of the transition? If the transition is of sec- 
ond order, what are the exponents? Are the exponents topology dependent? Is the 
transition temperature topology dependent? 



N 

- fiA, ■ nA,) - a ■ E • (2) 



i=0 
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2 Simulation Method and Autocorrelations 

First let us study the effect of tlie extrinsic curvature on tlie auto-correlation time 
of several observables, i.e., we are interested in the dynamics of a random surface 
induced by a Monte Carlo process. Such a study is of course a prerequisite for a 
detailed numerical study of the apparent transition which such surfaces undergo 
as a function of a bending rigidity. From our results which we present below we 
must conclude that extensive computer resources must be applied because of the 
very long relaxation times which even for moderate surface sizes can reach several 
tenthousand sweeps. 

We want to look at the closed dynamically triangulated random surfaces with- 
out self- avoidance. As models for such surfaces we take the first two topologically 
closed surfaces: The sphere and the torus (c.f. Figure [|). The torus is specified by 
identifying edges of the parameter space V whereas the sphere, because of its genus 
needs a slightly more complicated set-up procedure. 



Figure 1: Shown are two examples of configurations of dynamically triangulated 
random surfaces. The left picture shows a sphere and the right part a torus. 

What we are interested in is to calculate the auto-correlation function of an 
observable O function of time 



n (t m- < ^(o)^(^) >N-<o>^^ ... 



and its dependence on the number of nodes N of the surface. From the auto- 
correlation function we are then able to extract the auto-correlation time to{N). 
This is done calculating the integrated auto-correlation function, i.e., 

1 r+oo 

Tinw = 7; dt poit,N) . (4) 



2 7- 



oo 



We also compared the integrated auto-correlation time with the exponential correla- 
tion time and the statistical inefficiency, but found no disagreement within the errors. 
If the successively generated configurations can be considered as non-interacting 
beads as is done in the Rouse theory p7| , we will find for example for the radius 



of gyration 

r7^2 oc (5) 

and in general we expect 

To oc . (6) 

To calculate observables such as the gaussian part Tig of the Hamiltonian or the 
radius of gyration Rgyr, we use the standard Monte Carlo algorithm [29-31]. Such 
an algorithm induces a stochastic dynamics and all our results apply only to such 
an induced dynamics. One Monte Carlo sweep is completed when each triangula- 
tion point was given the chance for a displacement from its previous position and 
the edges of the triangulation were given the chance to re-connected or fiip to an 
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Figure 2: This figure demonstrates the elementary moves which are made during 
one update of a triangulated random surfaces. 

orthogonal position (c.f. Figure ^. The flip operation implements the sum over all 
possible triangulation) . This corresponds to one new configuration. 

Table |l| summarizes our results for the exponent a for the sphere. The exponents 
where extracted from the data of system sizes from 36 up to 288. At this point we 
should mention that the expected transition is at Ac ~ 1.5. Out data was collected 
above and below the transition. 

Above and below the transition the radius of gyration shows the Rouse behaviour. 
The gaussian part of the Hamiltonian shows clearly a dependence on the bending 
rigidity. Above the transition point the auto-correlation exponent changes to one 
and below it is almost zero. The same holds for the edge part of the Hamiltonian. 



A 


correlation time r of ... 


acceptance rate 


' '-gausstan 


' ^edqe 


flV" 


shift 


fiips 


0.00 


0.2 ±0.3 


0.2 ±0.2 


1.1 ±0.2 


0.51 


0.76 


1.00 


0.1 ±0.3 


0.4 ±0.2 


1.1 ±0.4 


0.51 


0.36 


1.75 


1.0 ±0.1 


0.8 ±0.2 


1.1 ±0.2 


0.47 


0.28-0.34 


3.00 


1.0 ±0.1 


1.0 ±0.1 


1.0 ±0.3 


0.50 


0.23-0.33 



Table 1: Table of the results for the exponent a on the auto- correlation times of 
dynamically triangulated random sphere. A is the measure of the bending rigidity. 
Shift and fiips give the acceptance probability for a displacement of a node and flip 
for an edge change. 

Almost all of the observations for the spherical case carry over to the case of the 
torus. 



3 Results of the Simulations 

Finite size scaling assumes, that there is only one relevant linear length scale, which 
is compared to the correlation length. To apply finite size scaling to the crumpling 
transition of dynamically triangulated random surfaces (DTRS) one must therefore 
assume a single length scale determined by the number of nodes and the internal 
dimension d of the surface 



Lo^N^''^ . (7) 

This internal dimension d also depends on the external properties of the surface and 
A 

So let us first look at the specific heat. If the transition is of second order we 
would have 
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where C* is a scahng function, which depends on how one implements the surface. 
At the critical A the scaling function is regular and the scaling hj^othese leads to 



N 



oc 



(9) 



for the scaling of the peak in the specific heat. If we assume a first order transition 
diverges as L'^ because of the ^-distribution of 



then C^""" 



34 



An evaluation of the specific heat C of DTRS (neglecting the metric contribution 
Hm) gives the following expression 



all 



D X' 



<n:> - <ne>' 



(10) 



The first part is related to the Gaussian Hamiltonian Tig and the second to the 
specific heat C of the edge extrinsic curvature He- The specific heat for the edge 
extrinsic curvature is shown in Figure ^ for the two topologies considered. The 
interpolation was done by the method of Ferrenberg and Swendsen ^ using 
histograms of He ■ 



Figure 3: Specific heat C (edge extrinsic curvature part) of the sphere (left part) 
and of the torus (right part). 



Figure 4: Maximum of the specific heat C]^"^ of the torus (O) and the sphere (□) 

Using the data of the specfic heat obtained by applying the extrapolation method, 
we can get a very accurate estimate of the positions of the maxima and of the heights. 
Figure ^ shows C]^"^ of the torus and the sphere. A change to L"^ behaviour is very 
unlikely and for that reason the data strongly suggest a continuous phase transition 
in agreement with previous work [8,14-17,38,39]. From the data shown in Figure ^ 
we can obtain the following upper boundaries of critical exponents 

Sphere:— < 0.00 ± 0.04 Torus:— < 0.06 ± 0.02 (11) 
ud ud 



Using this simple method, we cannot distinguish a diverging specific heat with 
very small but positiv a, a logarithmic divergence = and a power law cusp 
as < 0, where denotes the exponent of the singular part of the specific heat. We 
will use the abbreviation a instead of a^. 

Following Fisher these three cases may be distinguished with a fit of the 
form 
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C(AA) =A--(a\-"-1) +B 

a ^ ' 



AA 



121 



Figure ^ shows such fits for a = 0.1 (power law), a = —0.01 (near logarithmic) 
and a = —1.2,-2.0 (power law cusp). The fits clearly favour a power law cusp 
of the specific heat with a large negative value of a, although we were not able to 
estimate the exponent a precisely. 



Figure 5: Plot of the specific heat C against x = ^ i^AX " — 1 j with a = 0.1 (upper 
left), —0.01 (upper right), —1.2 (lower left) and a = —2.0 (lower right). 



Our results are in contrast to the estimates of Renken and Kogut a/ud = 
0.14(3) for the sphere. They assumed scaling above = 72, but the more accurate 
data in Figure H shows that this assumption is not true. The decreasing slope in 
Figure ^ is also present in their Figure 5 although this is slightly masked by 
larger error barrs and less data points. Catteral et all. |^8[ reported an estimate of 
a/i^d = 0.185(50), but they used a different discretization based on 0^ graphs with 
a flip acceptance rate of only 0(1%) at Ac ~ 1.5 and larger correlation times of the 
edge extrinsic curvature than in our discretization. 

Another possibility to determine the order of the transition is the cumulant Vn 
of the edge extrinsic curvature 



V 



TV 



1 (^e^)Ar 



3 {n^'f 



(13) 



N 



which behaves quite differently at temperature driven first- and second-order tran- 



sitions 33 



1. and 2.order: VnL^^ '^=" f T^T.Rxed (14) 

2.order: V;v|_ ^=°° | T = T,(iV) (15) 

1-order: V^U 1 - T = T,(iV) (16) 

and E_ are the energies of the system above and below the transition. For a 
very weak first-order transition {E^ ^ E_) we also have VNlmin ~ 2/3. 

We computed Vat defined by equation (O) using again the method of Ferren- 



berg and Swendsen 0. The resulting figures show the predicted single peak 
minima with wings following equation (p!^. The finite size dependence of the min- 
ima of VNlmin shown in Figure ^ distinctly favours the asymptotic behaviour in 
equation ([l^) and therefore a continuous phase transition or a very weak first order 
transition. 

In general, finite size scaling predicts also a shift AA = A^ — A^ of the effective 
transition 'temperature' A^ proportional to N~^{= L~'^) for a first- and proportional 
to N~^^'^^(= L'^^") for a second-order transition. Unfortunatly A^ of dynamically 
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Figure 6: Finite size dependence of the minimum VN\min of the reduced cumulant 
of the edge extrinsic curvature. (O) denotes data of the sphere, (□) those of the 
torus and the arrow the large N limit 2/3 of a continuous transition. 

triangulated random surfaces is not known. For that reason, we have to use a two- 
parameter fit with unknowns and ud. But we can improve the rehabihty of 
this fit, because we know more. First, the fit has to be a straight hne (neglecting 
corrections to scahng) and we have hm^r^oo ^A(A^) = 0. Second, the shift AA(A^) 
is different for the specific heat and the cumulant in general. Therefore the slope of 
the fitted lines will be different in general, but we still have limTv^oo AA(A^) = for 
both observables. 

Figure 7: Best fit of y = AA against x = N~^^'^'^' for the sphere (left part) and the 
torus (right part). (O) denotes the specific heat C data, (□) those of the cumulant 
V. 

Figure |^ shows the best fits for the sphere and the torus. The estimates of the 
parameters are 

A;!° = 1.51 ± 0.04 , z/rf = 3.2±0.5 (17) 

for the sphere and 

A~ = 1.47 ±0.02 , i^rf = 2.5±0.5 (18) 

for the torus. 

We turn now to the order parameter itself. Common practise is to take C, = R/L 
{R is the typical radius, L the linear size of the membrane) to be a suitable order 
parameter [41-47]. Initially it was also defined as Rg{L) = (L {L ^ oo), with 



the linear size L of the hexagon and the radius of gyration RgOcY^ij \ ~ 
A suitable choice for the order parameter therefore is 

C := R'/N (19) 

with 

1 / ^ ^ ^ 2\ 

CTj denotes the number of nearest neighbours, i.e. the number of links connected 
to a node i. The associated susceptibility is 

XR^ := L' {{C) - {(f) = ^ {{R') - {R'f) (21) 

Figure D shows the order parameter ( and Figure ^ the susceptibilty of the sphere. 
The results for the torus are similar. 
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Figure 8: Order parameter ( = /N of the sphere. Errors are smaller than symbol 



size. 



Figure 9: Susceptibility xr^ of the sphere (Lines to guide the eye). 
We also measured another possible order paramter C' = {R''^)/N, 

which exhibits a small increase near the phase transition and a slower decay of (' 
for A— >-0. The difference is caused by a change of the internal geometry near the 
phase transition [3^. 



With the data for the sphere in Figure |g and the corresponding data of the 
torus one can estimate the critical exponent (3/vd oi the order parameter assuming 
a second order transition. Here we use as the effective critical temperature the 
position of the peak of the specific heat. C^""^ (c.f. Figure p!0|). 

Figure 10: Scaling of the radius of gyration squared {R^{N)) at the position of the 
maximum of the specific heat (7™"^. (O) denots the data of the torus, (□) data of 
the sphere. 

An estimate of (3/vd using 

{R^N.X'Jf)) (xN^'^""^^ + ^ (22) 

results in 

Torus: [3/ud = 0.28 ± 0.02 , Sphere: p/ud = 0.35 ± 0.04 (23) 
and with the values of ud in equation (|l^) and (|18D 



Torus: p = 0.7 ± 0.2 , Sphere: ^ = 1.1 ± 0.2 (24) 

Figure |Tl] shows the scaling of the maxima of the susceptiblity x (equation (pTj) , 
Figure |^) . 
The results are 

Torus: = 0.62 ± 0.06 , Sphere: = 0.66 ± 0.06 (25) 

and, using the values of ud in equation (|1^ and (pTsP, 



Torus: 7 = 1.6 ± 0.5 , Sphere: 7 = 2.1 ± 0.5 . (26) 

With these large error bars and estimated values a ~ —1.5 the results are almost 
compatible with the scaling relation 
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Figure 11: Scaling of the susceptibility x (equation ^^), Fig. ^ in the case of a 
torus {O) and a sphere (□). 



a + 2/3 + 7 = 2 (27) 
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